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proposed  dredged  material  disposal  sites. 
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modeling  hurricane  storm  surge.  The 
hydrodynamic  model  used  was  the  ADCIRC- 
2DDT  code,  which  is  based  on  a  two-dimen¬ 
sional,  depth  integrated,  finite  element  for¬ 
mulation.  Hurricane  wind  stress  and  pressure 
forcing  from  Hurricane  Kate  (the  1985  histori¬ 
cal  storm  modeled)  were  produced  by  the 
HURWIN  code,  a  vertically  averaged, 
planetary  boundary  layer  wind  model.  The 
resulting  storm  surge  elevations  predicted  at  8 
stations  along  the  Florida  coast  correlated 
well  with  measured  storm  surge  elevations. 
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Conversion  Factors, 
Non-SI  to  SI  Units  of 
of  Measurement 


Non-SI  units  of  measurement  used  in  this  report  can  be  converted  to  SI 
units  as  follows: 


Multiply 

By 

To  Obtain 

knots 

0.5144444 

meters  per  second 

miles  per  hour 

1.609347 

kilometers  per  hour 

millibars 

1.0 

kilopascals 

Summary 


This  report  investigates  the  use  of  large  domains  in  modeling  hurricane 
storm  surge.  The  hydrodynamic  model  used  in  this  study  is  the 
ADCERC-2DDI  code  which  is  based  on  a  two-dimensional,  depth-integrated, 
finite  element  formulation.  Hurricane  wind-  stress  and  pressure  forcing  from 
Hurricane  Kate  are  produced  by  the  HURWIN  code,  a  vertically  averaged 
planetary  boundary  layer  wind  model.  Storm  surge  predictions  are  conducted 
over  three  computational  domains  each  having  a  different  size.  The  smallest 
domain  is  situated  on  the  continental  shelf,  another  domain  covers  the  Gulf  of 
Mexico,  and  the  largest  domain  includes  the  Western  North  Atlantic  Ocean, 
the  Caribbean  Sea,  and  the  Gulf  of  Mexico.  The  largest  domain  is  shown  to 
be  optimal  for  use  with  storm  surge  models.  The  influence  of  an  inverted 
barometer  condition  applied  at  the  open  boundary  is  examined  over  each 
computational  domain. 

Both  tidal  and  storm  singe  elevations  are  computed  over  the  large  scale 
domain.  Tidal  elevations  compare  well  with  measured  tidal  data  at  77  stations. 
Storm  surge  predictions,  which  include  both  tidal  and  inverted  barometer  open 
boundary  forcing,  show  excellent  correlation  with  measured  storm  surge  values 
when  compared  at  eight  stations  along  the  Florida  coast 
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1  Introduction 


The  hurricane  season  c  1992  included  such  destructive  and  well  publicized 
storms  as  Hurricane  Andrew  and  Hurricane  Iniki.  Apart  from  the  number  of 
lives  lost  during  these  hurricanes,  the  costs  due  to  material  damage  and 
disruption  to  local  economies  totalled  into  the  billions  of  dollars.  The 
devastation  from  these  recent  storms  reinforces  the  fact  that  coastal  flooding 
from  hurricane  storm  surge  remains  a  serious  problem  for  seaside 
communities.  In  the  interest  of  preventing  loss  of  human  life,  the  National 
Weather  Service  (NWS)  issues  warnings  when  hurricanes  threaten  U.S.  coastal 
areas.  NWS  hurricane  warnings  contain  information  concerning  the  range  of 
expected  peak  surge  heights  within  the  warning  area  (Jarvinen  and  Lawrence 
1985).  Computations  made  by  a  numerical  storm  surge  model  form  the  basis 
for  NWS  storm  surge  warnings  (Jelenianski  1979).  The  U.S.  Army  Corps  of 
Engineers  also  uses  storm  surge  models  to  design  coastal  protection  systems 
and  to  investigate  the  fate  of  dredged  material  disposed  in  the  coastal 
environment.  This  report  focuses  on  selection  of  an  appropriate  numerical 
storm  surge  model  and  its  use  in  simulating  hurricane  storm  surge. 

Storm  surge  is  a  long-period  wave  caused  by  extreme  wind  and  pressure 
forces.  Water  heights  associated  with  storm  surge  are  superimposed  on  water 
levels  generated  by  tidal  forcings.  The  period,  wavelength,  and  amplitude 
characteristics  of  a  storm  surge  depend  on  both  geometric  properties  of  the 
water  body  and  characteristics  of  the  meteorological  forcing  (Dendrou,  Moore, 
and  Myers  1985).  A  numerical  model  for  storm  surges  must  resolve  the 
physical  features  which  affect  storm  surge  generation  and  propagation.  Thus, 
an  appropriate  storm  surge  model  incorporates  complex  coastal  geometries, 
accounts  for  rapidly  varying  bathymetry  in  the  continental  slope  and  shelf 
regions,  and  permits  reasonable  boundary  condition  specification.  During  the 
course  of  development,  storm  surges  and  their  associated  forerunners 
encompass  a  large  range  of  spatial  scales.  An  acceptable  numerical  storm 
surge  model  must  represent  these  scales  of  motion  accurately  so  that  energy  in 
the  system  does  not  become  artificially  excited,  and/or  trapped  within  the 
system,  and/or  aliased  in  nonphysical  ways. 

A  finite  element  (FE)- based  numerical  model  with  its  inherent  grid 
flexibility  is  selected  as  the  ideal  formulation  for  a  storm  surge  model. 
Coastline  detail  is  easily  incorporated  and  nodal  densities  can  range  from  three 
to  four  orders  of  magnitude.  The  wide  variation  in  nodal  density  permits 
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significant  refinement  near  coastal  areas  and  in  regions  of  rapid  bathymetric 
changes  and  yet,  the  discrete  problem  remains  well  within  computational 
limits.  In  this  report,  an  FE-based  storm  surge  model,  ADCIRC-2DDI 
(Luettich,  Westerink,  and  Scheffner  1992),  which  Implements  the 
wave-continuity  equation  (WCE)  and  momentum  balance  equations,  is  applied. 
The  accuracy  of  these  equations  in  solving  various  shallow-water  problems  is 
well-documented  (Walters  1988;  Werner  and  Lynch  1989;  Walters  and  Wemer 
1989;  Gray  1989;  Foreman  1988;  Lynch  et  al.  1989;  Lynch  and  Werner  1991; 
Luettich,  Westerink,  and  Scheffner  1992). 

Aside  from  the  numerical  model,  careful  consideration  must  be  given  to  the 
specification  of  wind  stress  and  pressure  at  the  water  surface.  Wind  stress  and 
pressure  forcing  are  what  drive  the  numerical  storm  surge  model.  As  such, 
storm  surge  predictions  are  limited  by  the  accuracy  of  the  specified  wind  stress 
and  pressure  fields  (Hubbert,  Leslie,  and  Manton  1990;  Dendrou,  Moore,  and 
Myers  1985;  Flather  1984).  A  dynamic  numerical  model  of  the  planetary 
boundary  layer,  which  takes  advantage  of  available  meteorological  data  from 
historical  storms,  is  chosen  to  compute  the  wind  stress  and  pressure  forcing 
(Cardone,  Greenwood,  and  Greenwood  1992).  This  physically  based  wind 
model  offers  a  more  realistic  representation  of  the  meteorological  forcing  than 
more  common  empirical  forms  such  as  those  described  by  Cardone, 
Greenwood,  and  Greenwood  (1992)  and  those  used  by  Johns  et  al.  (1983a), 
Flather  (1984),  Jarvinen  and  Lawrence  (1985),  Hearn  and  Holloway  (1990), 
and  Westerink  et  al.  (1992). 

As  noted  by  Jarvinen  and  Lawrence  (1985),  the  initial  design  and  location 
of  the  basin  domain  is  a  critical  step  in  the  construction  of  a  storm  surge 
model.  Judicious  location  of  a  domain  results  in  optimal  storm  surge 
calculations.  Storm  surge  modeling  efforts  to  date  have  primarily  used 
domains  which  are  limited  to  the  continental  shelf  region  and  utilize  finite 
difference  methods  with  minimum  grid  flexibility  (Hearn  and  Holloway  1990; 
Al-Rabeh,  Eunay,  and  Cekirge  1990;  Hubbert,  Leslie,  and  Manton  1990; 
Lardner  and  Cekirge  1988;  Dube,  Sinah,  and  Roy  1986;  Dendrou,  Moore,  and 
Myers  1985;  Flather  1984;  Johns  et  al.  1983a, b).  Within  the  framework  of 
finite  differences,  the  need  for  grid  refinement  in  near  coastal  regions  is  often 
accommodated  by  a  nested  grid  approach.  A  small,  finely  discretized  grid  is 
contained  within  or  coupled  to  a  courser  regional  grid  (Hubbert,  Leslie,  and 
Manton  1990;  Al-Rabeh,  Eunay,  and  Cekirge  1990;  Dendrou,  Moore,  and 
Myers  1985;  Flather  1984).  Near  coastal  continuous  grid  compression  (Dube, 
Sinha,  and  Roy  1986;  Johns  et  al.  1983a, b)  is  another  technique  also 
employed.  While  this  technique  more  closely  approaches  the  grid  flexibility  of 
the  finite  element  method,  the  finite  element  method  remains  superior  in  its 
ability  to  accommodate  discretization  demands  in  an  optimal  fashioa 

In  this  report,  ADCIRC-2DDI,  a  WCE  FE-based  numerical  storm  surge 
model  is  used  to  investigate  the  role  of  the  computational  domain  in  accurate 
storm  surge  prediction.  Storm  surge  is  modeled  over  three  computational 
domains,  each  spanning  successively  larger  portions  of  the  study  area.  A 
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comparison  of  storm  surge  predictions  given  by  the  model  yields  two  major 
conclusions.  First,  a  computational  domain  encompassing  only  the  continental 
shelf  grossly  underpredicts  hurricane  storm  surge  and  generally  leads  to 
erroneous  values.  Secondly,  using  a  domain  which  extends  into  the  deep 
ocean  leads  to  the  most  accurate  representation  of  the  storm  surge.  On  the 
continental  shelf,  both  hurricane  winds  and  the  pressure  deficit  drive  surface 
water  up  against  the  shoreline.  Surface  water  elevations  in  the  deep  ocean, 
however,  are  mainly  affected  by  the  pressure  deficit.  Because  the  deep  ocean 
pressure  deficit  is  readily  obtained  and  since  the  storm  surge  model  is  driven 
by  surface  water  elevations  specified  at  the  open  ocean  boundary,  use  of 
surface  elevations  on  a  boundary  located  in  the  deep  ocean  yield  better  storm 
surge  predictions.  Tidal  elevations  at  the  open  boundary  in  a  large  domain 
model  can  also  be  easily  extracted  from  global  tidal  models.  The  interaction 
with  waters  adjacent  to  the  region  of  interest  is  also  important  when  modeling 
storm  surge.  A  large  domain  includes  these  adjacent  basins  allowing  all 
wavelengths  associated  with  a  storm  surge  to  propagate  throughout  the  domain 
and  onto  the  continental  shelf  where  development  of  storm  surge  is  most 
critical.  Thus,  a  large  domain  whose  open  boundary  lies  in  the  deep  ocean  is 
an  optimal  choice  for  storm  surge  prediction. 


Chapter  1  Introduction 


3 


2  Governing  Equations  and 
Numerical  Discretization 


The  computations  described  in  this  report  were  performed  using 
ADCIRC-2DDI,  die  depth-integrated  option  of  a  system  of  two-  and  three- 
dimensional  hydrodynamic  codes  named  ADCIRC  (Luettich,  Westerink,  and 
Scheffher  1992;  Westerink  et  aL  1992b).  ADCIRC-2DDI  uses  the  depth- 
integrated  equations  of  mass  and  momentum  conservation,  subject  to  the 
incompressibility,  Boussinesq,  and  hydrostatic  pressure  approximations.  Using 
the  standard  quadratic  parameterization  for  bottom  stress  and  neglecting 
baroclinic  terms  and  lateral  diffusion/dispersion  effects  leads  to  the  following 
set  of  conservation  statements  in  primitive  nonconservative  form  expressed  in  a 
spherical  coordinate  system  (Flather  1988;  Kolar  et  aL  1992): 


where 

t  =  time 

=  degrees  longitude  (east  of  Greenwich  positive)  and  degrees  latitude 
(north  of  the  equator  positive) 

=  free  surface  elevation  relative  to  the  geoid 
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U,V  =  depth-averaged  horizontal  velocities 
R  =  radius  of  the  Earth 
H  =  C  +  h  =  total  water  column 
h  =  bathymetric  depth  relative  to  the  geoid 
/  =  2Qsin$  =  Coriolis  parameter 
O  =  angular  speed  of  the  Earth 
p,  -  atmospheric  pressure  at  the  free  surface 
g  =  acceleration  due  to  gravity 
t|  =  effective  Newtonian  equilibrium  tide  potential 
p0  =  reference  density  of  water 
xa,  =  applied  free  surface  stress 

t.=  r  (V'  +  VT 
1  H 

Cf  =  bottom  friction  coefficient 

A  practical  expression  for  the  effective  Newtonian  equilibrium  tide 
potential  is  given  by  Reid  (1990)  as: 

11  (Mm)  -  E^VxW^)®08 


2*  ('Mo) 


LVA+v^) 


(4) 


where 

=  constant  characterizing  the  amplitude  of  tidal  constituent  of  species  j 
cty,  =  effective  Earth  elasticity  factor  for  tidal  constituent  of  species  j 
fM  =  time-dependent  nodal  factor 
vM  =  time-dependent  astronomical  argument 
j  =  0,1,2  =  tidal  species  (J  =  0,  declinational;  j  =  1,  diurnal;  j  =  2, 
semidiurnal) 

Lq  =  3siri24>  -  1 
L,  =  sin(2$) 

Li  =  cosfy) 

M  =  degrees  longitude  and  latitude,  respectively 
=  reference  time 

Tm  =  period  of  constituent  n  of  species  j 

Values  for  are  presented  by  Reid  (1990).  We  note  that  die  value  for  the 
effective  Earth  elasticity  factor  is  typically  taken  as  0.69  for  all  tidal 
constituents  (Schwiderski  1980;  Hendershott  1981)  although  its  value  has  been 
shown  to  be  slightly  constituent-dependent  (Wahr  1981;  Woodworth  1990). 

To  facilitate  an  FE  solution  to  Equations  1-3,  these  equations  are  mapped 
from  spherical  form  into  a  rectilinear  coordinate  system  using  a  Carte 
Parallelogrammatique  (CP)  projection  (Pearson  1990): 
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where 


x '  -  R(k-\)cos*0 

y'  *  *4> 


(5) 


(6) 


Xo4>0  =  center  point  of  the  projection 

Applying  the  CP  projection  to  Equations  1-3  gives  the  shallow  water  equations 
in  primitive  nonconservative  form  expressed  in  the  CP  coordinate  system: 

%  ,  cos4>0  *{UH)  ^  1  a(vffcos»)  _  n  (7) 

"57  cos$  dx'  cos<|>  a  y' 


at/  +  it  *u  +  vw_ 

~5F  cos$  3j^  dy/ 


u  +/ 


cos^0  a 


cos^  ax' 


Po 


Po* 


+  -T.U 


(8) 


av  +  cos<(>o  „  av  + 

"57  cos$  ax' 
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u+f 
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_a_  ^ 

3y'  Po 


+s(C 


+ 


(9) 


Utilizing  the  FE  method  to  resolve  the  spatial  dependence  in  the  shallow- 
water  equations  in  their  primitive  form  gives  inaccurate  solutions  with  severe 
artificial  near  2  x  Ax  modes  (Gray  1982).  However,  reformulating  tire 
primitive  equations  into  a  Generalized  Wave  Continuity  Equation  (GWCE) 
form  gives  highly  accurate,  noise-free,  FE-based  solutions  to  tire  shallow-water 
equations  (Lynch  and  Gray  1979;  Kinnmark  1984).  The  GWCE  is  derived  by 
combining  a  time-differentiated  form  of  the  primitive  continuity  equation  and  a 
spatially  differentiated  form  of  the  primitive  momentum  equations.  The  result 
is  recast  into  conservative  form,  reformulating  die  convective  terms  into 
nonconservative  form,  and  adding  the  primitive  form  of  the  continuity  equation 
multiplied  by  a  constant  in  time  and  space,  Xq  (Lynch  and  Gray  1979; 
Kinnmark  1984;  Luettich,  Westerink,  and  Scheffner  1992).  The  GWCE  in  the 
CP  coordinate  system  is: 
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32C  _  ac  .  a 

ap-  x°  "SF  a? -3F  [  "E^  a? 

-  v//i£  +  f*55l  v+/Vif  -  irf!5i±[i+f(c-ii) 

^  ^  *  yJ  cos*  ay[Po  sv* 


The  GWCE  (Equation  10)  is  solved  in  conjunction  with  the  primitive 
momentum  equations  in  nonconservative  form,  Equations  8  and  9. 

The  high  accuracy  of  GWCE-based  FE  solutions  is  a  result  of  their 
excellent  numerical  amplitude  and  phase  propagation  characteristics.  In  fact, 

Fourier  analysis  indicates  that  in  constant-depth  water,  by  using  linear 
interpolation,  a  linear  tidal  wave  resolved  with  25  nodes  per  wavelength  is 
more  than  adequately  resolved  over  the  range  of  Courant  numbers, 

C  -  \fghAt/Ax^l.O  (Luettich,  Westerink  and  Scheffher  1992).  Furthermore, 
the  monotonic  dispersion  behavior  of  GWCE-based  FE  solutions  avoids 
generating  artificial  near  2  x  Ax  modes,  which  plague  primitive-based  FE 
solutions  (Platzman  1981;  Foreman  1983).  The  monotonic  dispersion  behavior 
of  GWCE-based  FE  solutions  is  very  similar  to  that  associated  with  staggered 
finite  difference  solutions  to  the  primitive  shallow-water  equations  (Westerink 
and  Gray  1991).  GWCE-based  FE  solutions  to  the  shallow-water  equations 
allow  for  extremely  flexible  spatial  discretizations,  which  result  in  a  highly 
effective  minimization  of  the  discrete  size  of  any  problem  (Le  Provost  and 
Vincent  1986;  Foreman  1988;  Vincent  and  Le  Provost  1988;  Westerink  et  al.  1992). 
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The  details  of  ADCIRC,  an  application  of  the  GWCE-based  solution  to  the 
shallow-water  equations,  are  described  by  Luettich,  Westerink,  and  Scheffher 
(1992).  As  most  GWCE-based  FE  codes,  ADCIRC  applies  three-noded  linear 
triangles  for  surface  elevation,  velocity,  and  depth.  Furthermore,  decoupling  of 
the  time  and  space  discrete  form  of  the  GWCE  and  momentum  equations, 
time-independent  and/or  tri-diagonal  system  matrices,  and  Ml  vectorization  of 
all  major  loops  results  in  a  highly  efficient  code.  It  is  noted  that  wetting 
and/or  drying  of  elements  is  not  currently  accommodated  by  ADCIRC. 
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3  Specification  of  Hurricane 
Wind  Stress  and  Pressure 
Fields 


In  most  studies,  very  simple  descriptions  of  the  hurricane  wind  field  are 
used  to  drive  sophisticated  hydrodynamic  numerical  models.  This  disparity 
severely  limits  validation  of  these  coastal  ocean  models  and  puts  into  question 
the  usefulness  of  results  obtained  from  such  models.  In  order  to  reduce  the 
uncertainties  caused  by  a  poorly  specified  wind  and  pressure  field,  this  study 
incorporates  a  wind  model  based  on  the  basic  physical  and  dynamic 
characteristics  of  a  tropical  cyclone. 


Wind  stress  and  pressure  computations  throughout  this  report  were  carried 
out  by  a  modified  form  of  the  HURWIN  model.  The  HURWIN  model,  part  of 
the  Coastal  Modeling  System  at  the  U.S.  Army  Corps  of  Engineers  Waterways 
Experiment  Station,  is  based  on  the  vortex  model  of  Chow  (1971)  modified  by 
Cardone,  Greenwood,  and  Greenwood  (1992).  A  time  history  of  the  surface 
wind  field,  based  on  available  data  from  historical  hurricanes,  is  the  end 
product  of  the  HURWIN  model.  The  theoretical  basis  of  the  HURWIN  model 
is  the  equation  of  horizontal  motion,  vertically  averaged  through  the  depth  of 
the  planetary  boundary  layer  (PBL),  expressed  in  coordinates  fixed  to  the 
Earth: 


—  +fkxV  «  - 
dt 


_Lv/> 

P 


+  V-*tfVV  + 


\v\v 


(11) 


where 


d_ 

dt 


4-  +  V-V 

dt 


=  local  time  change  relative  to  fixed  coordinates 

V  =  two-dimensional  del  operator 

V  =  vertically  averaged  horizontal  wind  velocity 
/=  Coriolis  parameter 
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k  =  unit  vector  in  die  vertical  direction 
»  mean  air  density 

P  =  depth-averaged  pressure  in  the  PBL 
Kh  =  horizontal  eddy  viscosity  coefficient 
CD  s  surface  drag  coefficient 
h  =  depth  of  the  planetary  boundary  layer 

The  vertical  advection  of  momentum  is  assumed  to  be  small  compared  to  the 
horizontal  advection  and  is  thus  neglected.  Additionally,  die  shear  stress  at  the 
top  of  the  PBL  is  taken  to  be  zero. 

Total  pressure  is  prescribed  to  be  a  sum  of  the  pressure  associated  with  the 
tropical  cyclone  Pe  which  is  translating  with  the  storm  at  a  speed  Vt  and  a 

large-scale  background  pressure  P: 

P  -  Pt  +  P  02) 

The  background  or  far-field  pressure  can  be  specified  in  terms  of  a 
corresponding  geostrophic  flow  as: 


where  Vg  is  the  constant  geostrophic  velocity.  Substituting  in  the  pressure 
relationships  in  Equations  12  and  13  into  Equation  11,  the  integrated  wind 
field  within  the  PBL  becomes: 


(  \ 

( 1  'l 

.Lvp 

Ip  J 

Cn 

+  v-khvv  + 

D 

|V|V 


04) 


The  governing  equations  used  in  the  HURWIN  model  are  formed  by 
transforming  Equation  14  to  a  moving  Cartesian  coordinate  system  whose 
origin  is  located  at  the  low  pressure  center  of  the  storm. 

A  boundary  condition  appropriate  for  Equation  14  remains  to  be  defined. 
Noting  that  the  acceleration  and  horizontal  diffusion  of  momentum  at  the  edges 
of  the  computational  domain  are  neglected,  a  balance  is  implied  between  the 
Coriolis  force,  the  pressure  gradient  force,  and  the  surface  frictional  force  at 
these  boundaries. 

Specification  of  the  pressure  field  for  a  tropical  cyclone  Pc  is  accomplished 
using  the  well-known  exponential  pressure  law: 

P  -  P+hpe+»  05) 

*  *7# 
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where 


Ptyi  =  pressure  at  the  center  or  eye  of  the  storm 

Ap  -  pressure  anomaly,  P  - 

R  =  scale  radius,  often  assumed  equivalent  to  the  radius  to  the  maximum 
wind 

r  =  radial  distance  from  the  eye  of  the  stonn 

This  simple  but  widely  used  form  applies  to  a  radially  symmetric  pressure 
field.  An  option  is  available  in  which  the  pressure  anomaly  Ap  and  the  radius 
to  maximum  wind  R  can  be  specified  by  stonn  quadrant  With  this  option,  an 
asymmetrical  pressure  field  results  after  smoothing  between  the  pressure 
variations  specified  for  each  quadrant  The  theoretical  basis  of  die  PBL 
hurricane  wind  model,  HURWIN,  is  only  briefly  outlined  here.  A  more 
extensive  description  is  given  by  Cardone,  Greenwood,  and  Greenwood  (1992). 

The  computational  scheme  employed  by  the  HURWIN  model  involves  a 
nested  grid  system  consisting  of  five  rectangular  grids  each  with  a  uniform 
mesh  spacing.  This  mesh  spacing  doubles  for  successively  larger  grids.  For 
instance,  if  the  innermost  grid  has  a  spacing  of  5  km,  the  remaining  grids  have 
mesh  sizes  of  10  km,  20  km,  40  km,  and  80  km.  Thus,  the  maximum  areal 
extent  of  the  nested  grid  system  for  this  example  is  1,600  km2.  The  grid 
system  is  centered  at  the  eye  of  the  storm  and  moves  with  the  storm 
throughout  the  domain  as  dictated  by  the  form  of  the  governing  equations. 

The  numerical  discretization  of  the  governing  equations  is  based  on  a  finite 
difference  formulation,  the  details  of  which  are  given  by  Chow  (1971). 

The  HURWIN  model  is  based  on  the  concept  that  a  tropical  storm 
generally  changes  structure  relatively  slowly.  For  example,  the  structure  of  a 
hurricane  travelling  over  open  water  can  usually  be  well-represented  by 
parameters  specified  at  intervals  of  6  to  24  hr.  These  representative  states  are 
the  “snapshots"  of  the  wind  field.  Parameters  which  must  be  specified  for  a 
given  snapshot  include:  pressure  at  the  eye  of  the  storm,  far-field  pressure, 
radius  to  maximum  wind,  storm  track  direction,  speed  of  the  storm,  surface 
geostrophic  wind  velocity,  direction  of  the  surface  geostrophic  wind,  and  mesh 
spacing  of  the  innermost  grid.  Although  the  structure  typically  changes 
slowly,  storm  position  can  change  relatively  quickly.  Differences  in  the  time 
scale  between  storm  structure  and  position  changes  are  accommodated  by 
specifying  the  storm  center  at  every  hour.  The  hourly  information  needed 
includes  the  latitude  and  longitude  of  die  storm  center  and  the  relative 
weighting  between  snapshots  which  bracket  a  given  hour.  The  snapshot  and 
hourly  data  required  for  operation  of  the  HURWIN  model  are  obtained  from 
commonly  available  meteorological  data  for  historical  hurricanes.  The  full 
time  history  of  the  surface  wind  field  is  then  obtained  by  linear  interpolation  in 
time  using  this  series  of  characteristic  wind-field  snapshots.  The  snapshots  are 
calculated  at  discrete  times  during  the  course  of  the  storm  using  steady  state 
assumptions. 
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Despite  the  more  realistic  theoretical  basis  of  the  H  UR  WIN  model,  several 
limitations  arise  in  the  model's  ability  to  represent  the  hurricane  wind  field. 

The  most  prominent  shortcoming  relates  to  the  lack  of  dissipation  in  the  wind 
field  of  a  hurricane  making  landfall.  As  a  hurricane  approaches  land,  it 
generally  slows  and  winds  dissipate.  Consequently,  winds  on  the  right  side  of 
the  hurricane  are  somewhat  reduced  from  full  strength  as  they  encounter  the 
land.  As  the  winds  rotate  over  the  land,  their  speed  is  further  reduced  by 
frictional  forces  due  to  the  roughness  of  the  land.  Thus,  winds  on  the  left  side 
of  the  hurricane  coming  off  the  land  are  significantly  less  in  magnitude  than 
they  were  upon  initially  encountering  land.  The  PBL  wind  model  does  not 
represent  these  processes  associated  with  a  hurricane  making  landfall.  This 
limitation  of  the  PBL  hurricane  wind  model  leads  to  drying  of  coastal  elements 
on  the  left  side  of  the  hurricane  and  overprediction  of  peak  surge  on  the  right 
side  of  the  hurricane. 

Note  that  the  HURWIN  model  does  accommodate  terrains  of  varying 
roughness.  Upon  incorporation  of  surface  roughness  variations,  the  winds  over 
land  are  reduced  by  a  constant  factor.  But,  as  the  winds  come  off  the  land  and 
into  the  open  water,  the  wind  speed  will  uniformly  increase  in  accordance  with 
the  reduced  friction  of  the  open  ocean.  This  is  not  a  physically  realistic 
representation  of  the  wind  field  in  a  hurricane  as  it  approaches  landfall.  Left- 
and  right-ride  differences  in  the  wind  structure  are  not  addressed  by  this 
feature  in  the  PBL  hurricane  wind  model.  For  simulations  presented  in  this 
report,  all  surface  roughness  parameters  pertain  to  the  open  ocean. 

One  final  option  available  in  the  HURWIN  model,  which  is  not  utilized,  is 
the  specification  of  an  asymmetrical  storm  pressure  field.  Hurricanes  in  this 
report  are  represented  by  radially  symmetric  pressure  profiles.  Lack  of 
detailed  data  precludes  consideration  of  an  asymmetrical  hurricane  pressure 
representation. 

As  mentioned  previously,  the  HURWIN  model  has  been  slightly  modified 
so  as  to  directly  interface  with  the  ADCIRC-2DDI  hydrodynamic  model.  One 
addition  is  that  wind  speeds  computed  over  the  nested  grid  system  are 
converted  to  surface  wind  stresses  using  the  relationship  proposed  by  Garratt 
(1977): 

Il  =  CD£±|W|W  (16) 

Po  Po 

and 

ll  =  CD*±\W\Wy  (17) 
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where 


=  wind  stress  in  the  x  and  y  directions,  respectively 
PdPw  ~  ratio  of  air  density  to  average  density  of  seawater,  0.001293 
CD  -  frictional  drag  coefficient,  computed  (0.75  +  0.067 WOO. 001 
|W|  =  magnitude  of  wind  velocity 

Wx  ,Wy  =  components  of  the  wind  velocity  vector  in  the  x  and  y  directions, 
respectively 

In  die  H  UR  WIN  model,  pressure  is  not  explicitly  computed,  but  rather  the 
pressure  gradient  is  determined.  The  ADCIRC  model,  however,  requires  as 
input  a  pressure  field  expressed  as  an  equivalent  height  of  water 
Thus,  the  HURWIN  model  has  been  adjusted  to  produce  a  field  of  /Y(p»g) 
over  the  nested  grids.  Both  wind  stresses  and  P/(pj!)  values  are  linearly 
interpolated  directly  from  the  PBL  nested  grids  onto  the  finite  element 
computational  grid  used  by  the  ADCIRC  model  All  modifications  to  the 
HURWIN  model  result  in  an  output  of  hourly  wind  stress  and  P/ipJt)  values 
at  all  nodal  points  in  the  ADCIRC  computational  grid. 
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4  Hurricane  Kate  Description 


For  the  storm  surge  domain  studies  conducted  in  this  report.  Hurricane 
Kate  is  the  historical  storm  modeled.  The  track  of  Hurricane  Kate  through  the 
Western  Atlantic  Ocean  and  the  Gulf  of  Mexico  is  shown  in  Figure  1  for  1-hr 
increments.  Beginning  in  the  deep  Western  Atlantic  Ocean,  Hurricane  Kate 
was  recognized  as  a  tropical  storm  at  1800  GMT  on  November  15,  1985.  The 
storm  was  upgraded  24  hr  later  to  a  hurricane  as  winds  reached  over  80  mph.1 
A  minimum  pressure  of  953  mb  and  a  maximum  sustained  wind  speed  of 
about  121  mph  were  recorded  26.5  hr  prior  to  landfall  as  Kate  moved 
northeastward  through  the  Gulf  of  Mexico.  Landfall  of  Hurricane  Kate 
occurred  at  2230  GMT  on  November  21,  1985,  near  Panama  City,  Florida. 

The  tracking  of  Hurricane  Kate  continued  through  1800  GMT,  November  23, 
1985,  until  the  storm  had  significantly  dissipated  and  was  classified  as  an 
extratropical  storm.  Extensive  meteorological  analysis  of  Hurricane  Kate  as 
well  as  hydrographs  depicting  the  storm  surge  analysis  at  ten  stations  along  the 
Florida  coast  in  the  vicinity  of  Panama  City  can  be  found  in  the  report  by 
Garcia  and  Hegge  (1987). 

Wind  stress  and  pressure  forcing  for  the  ADCIRC-2D  simulations  ate 
computed  using  the  HURWIN  model.  Beginning  at  1800  GMT,  November  15, 
1985,  and  ending  8  days  later  at  1800  GMT,  November  23, 1985,  193  total 
hours  of  Hurricane  Kate  are  simulated.  During  this  period,  the  forward  speed 
of  Hurricane  Kate  ranged  from  3  to  26  knots.  The  far  field  pressure  of 
1,013  mb  was  taken  to  be  constant  throughout  the  simulation.  The  radius  to 
maximum  wind  varied  from  a  maximum  of  41  km  during  hurricane 
development  to  a  low  of  9  km  during  peak  intensity  of  the  storm.  Finally,  the 
innermost  nested  grid  spacing  is  set  at  5  km,  which  leads  to  1,600  km2  of  areal 
coverage  by  the  nested  grids  about  the  eye  of  the  hurricane. 


1  A  table  of  factors  for  converting  non-SI  units  of  measurement  to  SI  units  is  presented  an  page 
vii. 
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Figure  1 .  Track  of  Hurricane  Kate  through  the  Western  North  Atlantic  Ocean  and  into  the  Gulf 
of  Mexico  from  November  15, 1985,  1800  GMT  to  November  23, 1985, 1800  GMT 
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5  Comparison  of  Storm  Surge 
Predictions  Using  Three 
Domains 


A  numerical  storm  surge  model  is  comprised  of  several  major  components: 
the  governing  equations,  a  numerical  discretization,  wind  and  pressure  forcing, 
and  a  computational  domain.  All  of  these  components  can  significantly  affect 
the  ability  of  a  model  to  make  accurate  predictions.  Westerink  et  al.  (1992) 
previously  demonstrated  the  capabilities  of  the  ADCIRC-2D  model  with 
regards  to  the  governing  equations  and  numerical  discretization.  In  this  study, 
the  role  of  the  computational  domain  in  storm  surge  prediction  is  investigated 
through  consideration  of  three  different  domains.  Each  domain  telescopes 
outward  from  the  landfall  location  of  Hurricane  Kate  and  covers  successively 
larger  areas. 


Modeling  Strategy 

As  mentioned  previously,  most  storm  surge  modeling  efforts  to  date  have 
used  domains  which  extend  only  over  the  continental  shelf  region.  Physically, 
though,  appreciable  development  of  the  storm  surge  occurs  as  it  moves  up  the 
continental  slope  and  onto  the  shelf.  This  suggests  that  shelf  models  are  not 
optimal  for  storm  surge  predictioa  In  the  work  of  Westerink  et  al.  (1992),  a 
domain  comprised  of  the  entire  Gulf  of  Mexico  is  used.  Certainly  this  domain 
is  an  improvement  over  one  restricted  to  the  continental  shelf,  but  the  Gulf  of 
Mexico  is  plagued  by  resonant  frequencies  well  documented  by  Bungapong, 
Reid,  and  Whitaker  (1985).  These  resonant  modes  are  easily  excited  by  the 
shape  and  size  of  the  domain,  artifacts  in  the  domain  discretization,  forcing 
functions,  and/or,  the  boundary  conditions  and  can  interfere  with  an  accurate 
representation  of  a  hurricane  storm  surge.  A  very  large  domain  which  extends 
into  the  deep  ocean  is  ideal.  A  hurricane  can  progress  through  the  domain 
generating  and  propagating  storm  surge  in  a  more  natural  and  realistic  fashion. 
An  extensive  domain  more  fully  incorporates  the  factors  affecting  storm  surge 
generation  (i.e.,  the  bathymetry,  wind  and  pressure  forcing,  and  die  influence 
of  adjacent  basins).  In  the  deep  ocean,  changes  in  surface  water  elevation 
relate  mainly  to  pressure  deficits  and  not  wind  forcing,  so  for  the  large  domain 
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model,  an  inverted  barometer  condition  applied  at  the  open  boundary  is 
reasonable.  An  inverted  barometer  P/(pjz)  is  simply  die  height  to  which  the 
seawater  will  rise  due  to  pressure  forcing.  In  shallow  waters,  storm  surge  is 
generated  by  a  complex  interaction  between  bathymetry,  coastline  geometry, 
and  wind  forcing.  Thus,  an  inverted  barometer  condition  applied  over  shallow 
waters  does  not  accurately  represent  surface  water  elevations  due  to  hurricane 
forcing.  Any  errors  due  to  the  erroneous  specification  of  boundary  forcing  due 
to  a  storm  are  minimized  for  the  large  domain  model  since  the  open  ocean 
boundary  is  located  quite  far  from  the  coastal  region  of  interest.  Finally,  for  a 
large  domain  model,  tidal  elevation  data  can  easily  be  obtained  from  global 
ocean  models. 


Computational  Domain  Descriptions 

The  largest  domain,  the  east  coast  domain  shown  in  Figure  2,  encompasses 
the  Western  North  Atlantic  Ocean,  the  Caribbean  Sea,  and  the  Gulf  of  Mexico. 
The  east  coast  domain  includes  an  open  ocean  boundary  which  extends  from 
Glace  Bay,  Nova  Scotia,  to  the  vicinity  of  Corocora  Island  in  eastern 
Venezuela  along  the  60°W  meridian.  All  other  boundaries  are  defined  by  the 
eastern  coastlines  of  North,  Central,  and  South  America.  Topography  within 
the  domain,  depicted  in  Figure  2,  includes  the  continental  shelf,  whose  depths 
range  from  an  imposed  minimum  of  3  m  to  7  m  (Westerink,  Luettich,  and 
Scheffner  1993)  to  130  m,  the  continental  slope  which  has  a  typical  depth 
range  of  130  m  to  3,000  m,  and  the  continental  rise  and  deep  ocean  where 
bathymetries  range  from  3,000  m  upwards  to  almost  8,000  m.  Bathymetry  in 
the  Gulf  of  Mexico  is  specified  to  be  in  accordance  with  the  detailed  database 
used  by  Westerink  et  al.  (1992). 

The  east  coast  domain  discretization,  shown  in  Figure  3,  is  based  on 
Grid  V3  created  by  Westerink,  Luettich,  and  Scheffner  (1993)  for  the  Western 
North  Atlantic  Tidal  model.  Changes  to  Grid  V3  include  refinement  by  a 
factor  of  two  in  the  interior  regions  of  the  Gulf  of  Mexico  and  the  deep  ocean. 
Also,  the  coastline  and  coastal  waters  surrounding  Panama  City,  Florida,  the 
region  where  Hurricane  Kate  made  landfall,  have  been  substantially  refined. 
The  approximate  nodal  spacing  for  the  east  coast  domain  discretization  varies 
according  to  region.  As  an  example,  the  grid  size  in  coastal  areas  is 
approximately  0.006°,  while  in  the  deep  ocean,  spacing  is  about  1.15°.  The 
east  coast  domain  discretization  is  constructed  with  22,711  nodes  and  41,709 
elements,  reasonable  numbers  considering  the  level  of  refinement  and  the  areal 
extent  of  the  domain.  Optimal  gridding  and  the  flexibility  of  the  FE  method 
make  such  a  large  computational  problem  manageable.  Resolution  of  the 
domain  can  be  expressed  using  an  Af2  wavelength-to-grid  size  ratio  A# 
which  is  a  measure  of  the  number  of  nodal  points  used  to  resolve  a  wavelength 
of  the  Af2  tide.  A  minimum  criterion  for  resolution  of  the  M2  tide  is  A, 
equal  to  12.5.  The  X,W2/A,  ratio  for  the  east  coast  domain  ranges  from 
approximately  15  to  1,336. 
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Figure  3.  East  coast  domain  discretization 
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The  second  computational  domain,  the  Gulf  of  Mexico,  is  shown  in 
Figure  4.  This  domain  is  extracted  directly  from  the  east  coast  domain  and 
includes  the  entire  Gulf  of  Mexico.  Within  the  Gulf  of  Mexico  domain,  two 
well-defined  open  ocean  boundaries  exist  One  is  constructed  across  the  Strait 
of  Florida  extending  approximately  from  Cape  Sable,  FL,  to  Havana,  Cuba. 

The  second  open  ocean  boundary  stretches  across  the  Yucatan  Channel  in  the 
vicinity  of  Cancun,  Mexico,  and  Cabo  San  Antonio,  Cuba.  The  bathymetry, 
displayed  in  Figure  4,  for  the  Gulf  of  Mexico  domain  is  exactly  the  same  as 
that  specified  for  die  Gulf  of  Mexico  region  within  the  east  coast  domain. 
Again  for  simulation  purposes,  a  minimum  depth  of  3  m  is  imposed  and  the 
maximum  depth  in  the  Gulf  of  Mexico  domain  is  just  under  3,800  m.  The 
Gulf  of  Mexico  domain  discretization  entails  6,325  nodes  and  11,441  elements 
and  is  presented  in  Figure  5.  Approximate  nodal  spacings  of  the  discretization 
range  from  a  maximum  of  0.64°  to  a  minimum  of  0.006°  and  the  resolution, 
given  by  values  of  varies  between  15  and  721  grid  points  per 
wavelength  of  the  M2  tide. 

The  final  domain  under  study,  the  Florida  coast,  pictured  in  Figure  6,  is  the 
smallest  and  is  confined  to  the  continental  shelf.  The  Florida  coast  domain  is 
cut  directly  out  of  the  Gulf  of  Mexico  domain  and  spans  the  coastline  from  St 
Marks,  FL,  to  the  Mississippi  Sound,  reaching  outwards  to  the  edge  of  the 
continental  shelf.  The  open-ocean  boundaries  can  be  delineated  by  three 
connected  segments.  Two  of  the  segments  project  into  the  open  water  from 
the  coastline,  a  western  boundary  segment  beginning  near  Biloxi,  MS,  and  an 
eastern  boundary  segment  intersecting  the  coastline  halfway  between  St  Marks 
and  Cedar  Key,  FL.  Both  of  these  straight  line  segments  connect  with  an 
arcing  open  ocean  boundary,  which  loosely  follows  the  edge  of  the  continental 
shelf,  approximately  100-200  m.  Bathymetry  values  given  in  Figure  6  range 
from  the  minimum  of  3  m  to  nearly  240  m.  The  Florida  coast  domain  has  a 
discretization  of  1,933  nodes  and  3,145  elements  (see  Figure  7)  and  a  'kMj/AI 
range  of  27  to  597  grid  points  per  wavelength. 

Properties  of  each  of  the  three  computational  domains  under  study  are 
summarized  in  Table  1.  The  region  of  interest,  when  considering  the  storm 
surge  from  Hurricane  Kate,  is  the  Florida  coast  surrounding  the  point  of 
landfall  near  Panama  City,  FL.  It  is  noted  that  the  NWS  uses  a  standard 
semicircular  domain,  the  Pensacola  Bay  basin  (Jarvinen  and  Lawrence  1985), 
which  corresponds  roughly  to  the  Florida  coast  domain  used  in  this  study. 
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Figure  4.  Bathymetry  contours  in  increments  of  25,  50, 100,  200,  500, 1,000,  2,000,  and 
3,000  m  over  the  Gulf  of  Mexico  domain 
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Figure  5.  Gulf  of  Mexico  domain  discretization 


22 


Chapter  5  Comparison  of  Storm  Surge  Predctions  Using  Three  Domains 


Figure  6.  Bathymetry  contours  in  increments  of  25,  50, 100,  and  200  m  over  the  Florida  coast 
domain 


Figure  7.  Florida  coast  domain  discretization 
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Storm  Surge  Simulations 

Storm  surge  simulations  are  conducted  over  all  three  computational 
domains  using  the  ADCIRC-2D  code.  A  series  of  two  simulations  are  run, 
one  using  wind  and  pressure  forcing  and  the  other  applying  wind  and  pressure 
forcing  as  well  as  an  inverted  barometer  condition  on  the  open  boundary.  In 
neither  of  these  two  sets  of  simulations  is  tidal  forcing  considered,  either  at  the 
boundaries  or  on  the  interior  of  the  domain.  Identical  model  parameters  for  all 
simulations  are  selected  so  that  comparison  between  the  results  for  different 
domains  is  facilitated.  The  Florida  coast  and  Gulf  of  Mexico  domain 
simulations  exhibit  instabilities  due  to  drying  of  elements,  a  feature  not 
addressed  by  the  ADORC  model.  Consequently,  the  model  is  run  using  only 
the  nonlinear  bottom  friction  option  (i.e.,  the  convective  and  finite  amplitude 
terms  are  not  included).  A  constant  value  of  0.003  is  specified  as  the  Manning 
friction  coefficient,  Cp  throughout  all  domains.  Horizontal  momentum 
diffusion  is  parameterized  using  an  eddy  viscosity  formulation  where  the  eddy 
viscosity  coefficient  is  set  equal  to  5  m2/sec.  A  GWCE  parameter  which 
represents  the  balance  between  the  primitive  continuity  and  wave  equation 
portion  of  the  GWCE,  is  defined  equal  to  0.001. 

Simulations  were  spun  up  from  homogeneous  initial  conditions  using  a 
1-day  ramp  in  time.  Application  of  the  hyperbolic  ramp  function  reduces  the 
excitation  of  nonphysical  short  wavelength  frequencies.  An  identical  ramp 
function  of  1  day  length  is  applied  to  the  wind  and  pressure  forcing  as  well  as 
to  an  inverted  barometer  boundary  condition  when  used.  Actual  simulations 
were  run  over  8.25  days,  a  period  which  includes  the  1-day  rampup.  During 
the  first  6  hr  of  the  simulation,  initial  hurricane  winds  and  pressure  forcing  are 
held  stationary.  Following  these  6  hr,  storm  surge  predictions  begin  at  1800 
GMT,  November  15,  1985,  and  continue  for  8  days  ending  1800  GMT, 
November  23,  1985.  A  time-step  of  45  sec  is  used  throughout  the  simulation 
period. 

Elevation  stations  have  been  placed  throughout  all  three  domains. 

Locations  of  the  elevation  stations  are  identical  in  overlapping  regions  between 
the  domains.  Eleven  stations  are  placed  along  the  Florida  Coast  between  St 
Marks,  FL,  and  the  Mississippi  Sound,  and  fifteen  stations  are  located  along 
the  open  ocean  boundary  of  the  Florida  Coast  domain.  Around  the  remaining 
coastline  of  the  Gulf  of  Mexico  domain  there  are  16  elevation  stations,  and  12 
more  stations  are  placed  along  the  open  ocean  boundary  of  die  Gulf  of  Mexico 
domain.  Finally,  1 1  additional  elevation  stations  are  specified  in  the  east  coast 
domain,  one  at  Havana,  Cuba,  and  the  remaining  10  on  the  open  ocean 
boundary.  All  65  elevation  station  names  and  locations  are  shown  in  Figures 
8,  9,  and  10  for  each  computational  domain. 
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Figure  9.  Names  and  locations  of  elevation  stations  in  the  Gulf  of  Mexico  domain 
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Two  series  of  simulations  are  conducted,  each  implementing  a  different 
open  boundary  condition.  In  the  first  series  of  simulations,  zero  elevation  is 
specified  at  all  open  boundary  nodes  in  each  domain.  The  second  series  of 
simulations  applies  an  inverted  barometer  condition  at  all  open  boundary 
nodes.  For  these  simulations,  no  calibration  or  tuning  of  pararr  ters  was 
performed,  either  in  the  weather  model  or  in  the  hydrodynamic  model. 

Several  observations  arise  from  a  comparison  of  storm  surge  predictions 
over  each  domain.  The  storm  surge  elevations  predicted  at  St.  Marks,  Turkey 
Point,  and  Apalachicola,  FL,  are  shown  in  Figures  11a,  12a,  and  13a  using  a 
zero  elevation  open  boundary  condition  and  in  Figures  lib,  12b,  and  13b 
using  an  inverted  barometer  open  boundary  condition.  Storm  surge  profiles  at 
these  stations  are  representative  of  conditions  on  the  right-hand  side  of 
Hurricane  Kate  along  the  Florida  coast.  Gearly,  the  results  point  out  the  fact 
that  the  Florida  coast  domain  significantly  underpredicts  peak  storm  surge  and 
generally  models  the  storm  surge  rather  poorly  as  compared  to  the  other 
domain  results.  For  elevation  stations  near  Panama  City,  Destin,  and 
Pensacola,  FL,  stoim  surge  profiles  on  the  left-hand  side  of  Hurricane  Kate  are 
presented  in  Figures  14a,15a,  and  16a  for  the  zero  elevation  boundary 
condition  and  in  Figures  14b,  15b,  and  16b  for  the  inverted  barometer 
condition.  When  the  hurricane  winds  come  off  the  land,  the  Florida  coast 
domain  simulations  exhibit  excessive  negative  surges  indicating  that  drying  of 
coastal  elements  is  occurring.  Storm  surge  predictions  from  these  simulations 
lead  one  to  conclude  that  the  Florida  coast  domain  is  an  inadequate  domain  for 
simulating  storm  surge  when  compared  with  results  from  the  two  other 
domains  considered  here.  Recall  that  the  Florida  coast  and  the  Gulf  of  Mexico 
domains  were  cut  from  the  east  coast  domain  such  that  the  discretization  for 
all  domains  remain  the  same.  Thus,  any  differences  in  storm  surge  prediction 
over  the  various  domains  are  due  strictly  to  the  size  of  the  domain  and  the 
boundary  conditions  used  to  drive  the  flow.  Reiterating,  a  computational 
domain  which  covers  only  the  continental  shelf  is  inappropriate  for  storm  surge 
prediction.  This  is  related  to  the  significant  development  of  the  storm  surge 
which  occurs  as  a  hurricane  moves  over  the  rapid  bathymetric  changes  of  the 
continental  slope  and  onto  the  contineiual  shelf.  Currently,  most  storm  surge 
models  still  use  continental  shelf  domains,  clearly  an  inappropriate  domain. 
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Figure  1 1 .  Computed  storm  surge  tor  Hurricane  Kate  at  St.  Marks,  Florida 
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b.  Inverted  barometer  open  boundary  condition 
Figure  12.  Computed  storm  surge  for  Hurricane  Kate  at  Turkey  Point,  Florida 
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a.  Zero  elevation  open  boundary  condition 


b.  Inverted  barometer  open  boundary  condition 

Figure  13.  Computed  storm  surge  for  Hurricane  Kate  at  Apalachicola,  Florida 
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b.  Inverted  barometer  open  boundary  condition 


Figure  14.  Computed  storm  surge  for  Hurricane  Kate  at  Panama  City,  Florida 
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b.  Inverted  barometer  open  boundary  condition 


Figure  16.  Computed  storm  surge  for  Hurricane  Kate  at  Pensacola,  Florida 
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Storm  surge  predictions  for  the  Gulf  of  Mexico  and  east  coast  domains 
have  peak  surges  which  coincide  in  time  and  magnitude.  However, 
oscillations  plague  the  storm  surge  solution  over  the  Gulf  of  Mexico  domain  as 
clearly  indicated  in  all  figures  depicting  storm  surge  associated  with  the  zero 
elevation  boundary  condition.  These  oscillations  are  characteristic  of  the  Gulf 
of  Mexico  whose  resonant  modes  can  be  excited  by  the  forcing  functions, 
boundary  conditions,  and/or  a  specific  domain  as  well  as  the  domain 
discretization.  When  the  inverted  barometer  boundary  condition  is  used,  the 
oscillations  become  less  pronounced  but  still  are  evident.  While  in  general  the 
inverted  barometer  condition  leads  to  an  improvement  in  the  storm  surge 
predictions,  an  inverted  barometer  condition  is  not  all  that  is  needed  for 
accurate  storm  surge  prediction.  Figures  17-19  compare  inverted  barometer 
values  on  the  open  ocean  boundary  of  the  Florida  coast  domain  to  storm  surge 
elevations  generated  over  the  Gulf  of  Mexico  and  east  coast  domains.  From 
these  comparisons,  it  is  obvious  that  an  inverted  barometer  condition  near  the 
continental  shelf  break  will  grossly  underestimate  storm  surge  entering  the 
domain.  This  is  expected  since  the  storm  surge  in  shallow  water  is  created  by 
the  interaction  of  both  wind  and  pressure  forcing  with  local  bathymetry  and 
coastal  geometry.  Over  the  east  coast  domain,  the  inverted  barometer 
condition  leads  to  minimal  changes  in  storm  surge  elevation  as  shown  by 
Figure  20,  in  which  storm  surge  profiles  for  the  zero  elevation  and  inverted 
barometer  boundar  y  conditions  are  compared  at  St.  Marks,  FL.  In  the  large 
domain  model,  changes  of  surface  water  elevation  in  the  deep  ocean  are 
largely  due  to  pressure  deficits.  These  pressure  deficits  can  be  quite  small  in 
the  deep  ocean  and  thus  surface  water  elevations  specified  at  this  deep  ocean 
boundary  have  a  minimal  influence  on  storm  surge  generation  in  the  near 
coastal  region.  The  storm  surge  predictions  shown  in  Figure  20  are  evidence 
of  the  reduced  effect  of  the  deep  ocean  boundary  conditions. 

Storm  surge  predictions  over  three  different-sized  domains  were  compared 
and  the  largest,  the  east  coast  domain,  was  selected  as  the  optimal 
computational  domain.  A  large  domain  permits  a  realistic  response  of  the 
coastal  region  to  wind  and  pressure  forcing.  The  generation  and  propagation 
of  storm  surge  in  a  large  domain  is  well  represented  because  bathymetric 
changes  in  the  continental  slope  and  shelf  region,  complex  coastal  geometries, 
and  the  influence  of  adjacent  basins  are  properly  incorporated.  In  addition,  a 
large  domain  has  open  boundaries  in  the  deep  ocean  which  are  free  from  the 
complex  interactions  between  coastal  geometry,  bathymetry,  and  wind  and 
pressure  forcing  near  or  on  the  shelf.  Specification  of  open  boundary 
elevations  in  the  deep  ocean  are  thus  simplified.  Also,  since  the  deep  ocean 
boundary  is  far  from  the  coastal  region,  elevations  specified  at  this  boundary 
have  a  minimal  influence  on  coastal  storm  surge  generation. 
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ure  17.  Comparison  of  computed  storm  surge  for  Hurricane  Kate  using  an 
inverted  barometer  open  boundary  condition  at  Station  14  located 
on  the  Florida  coast  domain  open  boundary 
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Figure  18.  Comparison  of  computed  storm  surge  for  Hurricane  Kate  using  an 
inverted  barometer  open  boundary  condition  at  Station  20  located 
on  the  Florida  coast  domain  open  boundary 
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19.  Comparison  of  computed  storm  surge  for  Hurricane  Kate  using  an 
inverted  barometer  open  boundary  condition  at  Station  23  located 
on  the  Florida  coast  domain  open  boundary 


Figure  20.  Comparison  of  computed  storm  surge  for  Hurricane  Kate  over  the 
east  coast  domain  using  both  a  zero  elevation  and  an  inverted 
barometer  open  boundary  condition  at  St.  Marks,  Florida 
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6  Comparison  Between 
Computed  and  Observed 
Storm  Surge  Elevations 
Using  a  Large  Domain 


As  a  demonstration  of  the  storm  surge  prediction  capabilities  of  the 
ADCIRC  code.  Hurricane  Kate  is  modeled.  Computed  storm  surge  profiles  at 
eight  elevation  stations  are  compared  to  measured  storm  surge  data  for 
verification  of  the  model  results.  From  the  computational  domain 
comparisons,  the  east  coast  domain  is  selected  as  being  the  appropriate 
computational  domain  for  the  storm  surge  simulations. 

The  ADCIRC-2D  model  is  run  with  all  nonlinear  options  with  fire 
exception  of  the  finite  amplitude  terms.  These  terms  were  not  included  in  the 
formulation  because  of  problems  created  by  the  drying  of  coastal  areas  due  to 
the  excessive  wind  speeds  coming  off  the  land  on  die  left  side  of  the 
hurricane.  All  other  model  parameters,  including  the  time-step,  remain  as 
specified  for  the  domain  comparison  simulations  detailed  in  the  previous 
section. 

Another  addition  to  the  runs  presented  in  the  previous  section  is  that  tidal 
forcing  is  now  included  in  the  storm  surge  simulations.  The  tidal  potential 
function  is  applied  throughout  die  interior  of  die  domain  and  tidal  elevation 
forcing  is  specified  along  the  open  ocean  boundary.  Five  tidal  constituents,  Ku 
Ox,  P j,  M2,  and  S2,  are  considered  both  on  the  interior  and  at  the  open  ocean 
boundary.  Table  2  details  the  frequencies,  nodal  factors,  equilibrium 
arguments,  and  Earth  elasticity  factors  used  for  each  constituent  The  values 
of  these  tidal  parameters  are  defined  with  reference  to  November  19, 198S, 
0:00  GMT.  Along  die  open  ocean  boundary,  tidal  elevations  are  obtained  from 
Schwiderski's  global  model  (Schwiderski  1979;1981a-g).  For  these 
simulations,  the  amplitude  of  die  M2  tide,  as  reported  by  Schwiderski,  is 
reduced  by  10  percent  a  correction  in  accordance  with  the  error  analysis  of 
Westerink,  Luettich,  and  Scheffiier  (1993)  as  determined  by  their  tidal  model 
of  the  Western  Atlantic  Ocean. 
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Table  2 

Tidal  Potential  Conatanta  Uaed  for  Tidal  and  Storm  Surge 
Slmulatlona 

Constituent 

Period  (hr) 

Potential 

Amplitude 

(m) 

Earth  Tide 
>**1uctk>n 
i’i^tor 

Nodal 

Factor 

EquNferium  | 
Argument  1 

(degrees)  | 

K,  lun  t-solar 

23.934470 

0.141565 

0.736 

1.086 

321.3  | 

O,  principal  lunar 

25.819342 

0.100514 

0.695 

1.140 

271.35  1 

P,  principal  solar 

24.065890 

0.046843 

0.706 

1.000 

33.00  j 

M,  principal  lunar 

12.420601 

0.242334 

0.693 

0.974 

230.34  | 

S,  principal  solar 

12.000000 

0.112841 

0.693 

1.000 

0.00  ] 

Before  proceeding  to  storm  surge  simulation,  verification  of  the  tidal 
responses  throughout  the  east  coast  domain  is  undertaken.  A  45-day 
simulation  is  conducted  with  a  ramp-up  period  of  12  days.  Results  are 
recorded  alter  30  days  of  computation  to  ensure  dissipation  of  free  modes 
within  the  Gulf  of  Mexico  and  Caribbean  as  suggested  by  Westerink.  Luettich, 
and  Scheffner  (1993).  Computed  tidal  elevations  at  77  stations,  placed 
throughout  the  domain,  are  compared  to  measured  data  obtained  at  these  same 
stations.  Details  regarding  the  location  of  the  stations  and  the  source  of  the 
tidal  data  can  be  found  in  Westerink,  Luettich,  and  Scheffner  (1993).  Tidal 
predictions  agree  well  with  the  measured  data  as  seen  in  Figures  21-26.  The 
slight  downward  adjustment  of  the  M2  amplitude  specified  on  the  open 
boundary  produces  significantly  improved  tidal  elevations  at  all  stations  as 
compared  to  those  computed  by  Westerink,  Luettich,  and  Scheffiier  (1993). 

Storm  surge  predictions  over  the  east  coast  domain  include  the  tidal  forcing 
just  described  as  well  as  an  inverted  barometer  condition  on  the  open  ocean 
boundary.  The  simulation  begins  with  36.75  days  of  tidal  forcing  before  the 
hurricane  enters  the  domain.  Initial  hurricane  wind  stress  and  pressure  forcing 
are  then  held  stationary  for  6  hr.  Following  this  period,  storm  surge 
computations  continue  for  8  days  beginning  at  1800  GMT,  November  15, 

1985,  and  ending  1800  GMT,  November  23,  1985.  Storm  surge  elevations  are 
recorded  every  3.75  min  at  eight  elevation  stations  located  along  the  Florida 
coastline  from  St  Maries,  FL,  to  the  Mississippi  Sound.  Measured  storm  surge 
values  are  available  at  these  same  locations.  Comparisons  of  numerical  storm 
surge  predictions  with  measured  storm  surge  elevations  are  presented  in 
Figures  27-34.  These  results  indicate  that  the  storm  surge  from  Hurricane  Kate 
is  well  represented  by  the  ADCIRC-2D  model.  Overprediction  of  the  peak 
surge  and  the  significant  drop  in  surge  at  stations  to  the  west  of  Apalachicola, 
FL,  are  both  artifacts  of  the  HURWIN  model  used  to  generate  the  wind  stress 
forcing. 
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A  severe  limitation  of  the  HURWIN  wind  model  is  its  inability  to  account 
for  frictional  dissipation  of  the  winds  as  a  hurricane  makes  landfall  As  winds 
rotate  over  land  their  magnitude  remains  high  and  subsequent  drying  of  coastal 
regions  occurs.  This  drying  causes  a  large  reduction  in  storm  surge  elevation, 
so  much  so,  in  fact,  that  significant  negative  storm  surge  values  are  recorded 
which  are  not  in  accordance  with  measured  data.  The  HURWIN  hurricane 
wind  model  does  not  account  for  dissipation  of  hurricane  winds  either  in  the 
open  water  prior  to  landfall  or  as  winds  rotate  over  land  and  come  off  the 
shore.  Thus,  the  wind  stress  values,  computed  by  the  HURWIN  model  and 
used  by  the  ADC  IRC  model,  are  generally  too  high  in  the  region  of  hurricane 
landfall. 


Figure  21 .  Comparison  of  computed  and  measured  tidal  elevations  over  the 
east  coast  domain  at  Key  West,  Florida 
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Figure  22.  Comparison  of  computed  and  measured  tidal  elevations  over  the 
east  coast  domain  at  Cedar  Key,  Florida 


Figure  23.  Comparison  of  computed  and  measured  tidal  elevations  over  the 
east  coast  domain  at  St.  Marks  Light,  Florida 
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ure  24.  Comparison  of  computed  and  measured  tidal  elevations  over  the 
east  coast  domain  at  Alligator  Bayou,  Florida 


Figure  25.  Comparison  of  computed  and  measured  tidal  elevations  over  the 
east  coast  domain  at  Southwest  Pass,  Louisiana 
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Figure  26.  Comparison  of  computed  and  measured  tidal  elevations  over  the 
east  coast  domain  at  Western  Florida  Outer  Shelf 


Figure  27.  Comparison  of  computed  and  measured  storm  surge  elevations 
for  Hurricane  Kate  over  the  east  coast  domain  at  St.  Marks, 
Florida 
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Lire  28.  Comparison  of  computed  and  measured  storm  surge  elevations 
for  Hurricane  Kate  over  the  east  coast  domain  at  Turkey  Point, 
Florida 


NOVEMBER,  1985 


Figure  29.  Comparison  of  computed  and  measured  storm  surge  elevations 
for  Hurricane  Kate  over  the  east  coast  domain  at  Shell  Point, 
Florida 
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Figure  30.  Comparison  of  computed  and  measured  storm  surge  elevations 
for  Hurricane  Kate  over  the  east  coast  domain  at  Carrabelle, 
Florida 


Figure  31 .  Comparison  of  computed  and  measured  storm  surge  elevations 
for  Hurricane  Kate  over  the  east  coast  domain  at  Apalachicola, 
Florida 
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Figure  34.  Comparison  of  computed  and  measured  storm  surge  elevations 
for  Hurricane  Kate  over  the  east  coast  domain  at  Pensacola, 
Florida 
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7  Conclusions 


Most  storm  surge  modeling  efforts  to  date  have  used  domains  that  are 
confined  to  the  continental  shelf.  Little  effort  has  been  given  to  understanding 
the  effect  of  domain  size  on  storm  surge  predictions  and  the  subsequent 
selection  of  an  appropriate  storm  surge  domain.  In  this  report,  storm  surge 
calculations  over  three  successively  larger  domains  clearly  illustrate  that  storm 
surge  models  restricted  to  the  continental  shelf  severely  underestimate  the  peak 
surge  and  do  not  reproduce  coastal  storm  surge  profiles  in  an  acceptable 
fashion.  Storm  surge  intensification  results  from  the  complex  interactions 
between  coastal  geometry,  bathymetry,  and  wind  and  pressure  forcing  as  water 
is  pushed  up  the  continental  slope  and  onto  die  shelf.  This  report  demonstrates 
that  these  processes  cannot  be  properly  accounted  for  at  the  open  boundaries  of 
a  shelf  model.  Thus,  a  storm  surge  model  using  a  continental  shelf  domain  is 
inappropriate. 

Another  size  domain  is  one  that  covers  the  entire  Gulf  of  Mexico. 

Solutions  obtained  for  storm  surge  elevation  over  die  Gulf  of  Mexico  domain 
exhibit  oscillatory  behavior.  Such  oscillations  are  due  to  the  existence  of  free 
modes  in  the  Gulf  of  Mexico  basin.  These  resonant  modes  in  the  Gulf  of 
Mexico  are  well  documented  and  can  be  easily  excited  by  a  number  of  factors: 
numerical  disc  ization,  forcing  functions,  and/or  boundary  conditions.  Thus, 
storm  surge  elevations  obtained  over  the  Gulf  of  Mexico  are  considered 
unreliable  due  to  the  nature  of  die  basin. 

A  very  large  domain,  which  includes  the  Western  North  Atlantic  Ocean, 
the  Gulf  of  Mexico,  and  the  Caribbean,  is  found  to  be  ideal.  A  hurricane  can 
progress  through  the  domain  generating  and  propagating  storm  surge  in  a 
natural  and  realistic  fashion.  The  inclusion  of  basins  adjacent  to  the  coastal 
region  allows  all  wavelengths  associated  with  storm  surge  to  propagate  through 
the  domain  and  onto  the  continental  shelf  where  development  of  storm  surge  is 
most  critical.  Factors  that  affect  storm  surge  generation  and  propagation,  such 
as  bathymetric  changes  in  the  continental  slope  and  shelf  region,  complex 
coastal  geometries,  and  wind  and  pressure  forcing,  are  all  properly 
incorporated  into  a  large  domain  model.  Another  advantage  is  that  the  open 
boundaries  of  a  large  domain  lie  in  the  deep  ocean  and  are  far  from  the 
intricate  processes  occurring  in  the  continental  shelf  region.  In  the  deep  ocean, 
changes  in  the  surface  water  elevation  are  largely  due  to  pressure  deficits 
whereas  on  the  continental  shelf  both  pressure  and  wind  forcing  are  also 
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important  An  inverted  barometer  condition  can  be  easily  and  appropriately 
specified  for  the  deep  ocean  boundary.  Additionally,  tidal  elevations  at  the 
open  boundary  in  a  large  domain  model  are  readily  extracted  from  global  tidal 
models.  Again,  any  errors  due  to  the  erroneous  specification  of  boundary 
forcing  are  reduced  in  a  large  domain  model  because  the  open  ocean  boundary 
is  located  quite  far  from  the  coastal  region  of  interest.  For  all  of  these  reasons 
a  large  domain  whose  open  boundary  lies  in  the  deep  ocean  is  optimal  for 
storm  surge  prediction. 

The  finite  element  formulation  used  in  the  ADCIRC  hydrodynamic  model 
facilitates  use  of  a  large  computational  domain.  Flexibility  of  the  finite 
element  method  leads  to  easy  incorporation  of  coastline  detail  and  nodal 
densities  which  can  range  from  three  to  four  orders  of  magnitude.  A  wide 
variation  in  nodal  density  arises  due  to  significant  refinement  in  coastal  areas 
and  in  regions  of  rapid  bathymetric  changes  where  storm  surge  generation  is 
important  and  subsequent  coarse  discretizations  in  the  deep  ocean  where 
processes  occur  more  gradually  and  are  of  less  interest.  The  efficiency  of  the 
finite  element  method  leads  to  a  discrete  problem,  associated  with  the  large 
domain,  which  remains  well  within  computational  limits. 

Despite  the  more  realistic  theoretical  basis  of  the  HURW1N  model,  the 
model  does  not  represent  the  processes  associated  with  a  hurricane  making 
landfall.  Frictional  dissipation  of  the  winds  as  a  hurricane  makes  landfall  is 
not  addressed  in  either  die  open  water  prior  to  landfall  or  for  winds  that  rotate 
over  land  and  come  off  the  shore.  The  artificially  high  magnitude  of  the  wind 
field  leads  to  overprediction  of  peak  surges  and  drying  in  coastal  regions, 
features  not  supported  by  the  measurement  data.  Storm  surge  predictions  over 
the  large  domain  can  be  improved  by  incorporating  into  the  wind  model  the 
dissipative  processes  involved  when  a  hurricane  makes  landfall. 
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